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ON THE DETERMINANT FORMULA IN THE INVERSE 

SCATTERING PROCEDURE WITH A PARTIALLY KNOWN 

STEPLIKE POTENTIAL 

ODILE BASTILLE, ALEXEI RYBKIN 



Abstract. We are concerned with the inverse scattering problem for the full 
r ~t • line Schrodinger operator —d%. +q(x) with a steplike potential q a priori known 

on R+ = (0,oo). Assuming <?|m + is known and short range, we show that the 
unknown part g|n_ of q can be recovered by 



q\ R _ Or) = -29* logdet (l + (1 + WT+)- 1 G a! ) , 



where MJ is the classical Marchenko operator associated to q\m + and G x is a 
*S^ , trace class integral Hankel operator. The kernel of Gx is explicitly constructed 

£)_l in term of the difference of two suitably defined reflection coefficients. Since 

q\m_ is not assumed to have any pattern of behavior at —00, defining and 
analyzing scattering quantities becomes a serious issue. Our analysis is based 
upon some subtle properties of the Titchmarsh-Weyl m-function associated 
with R_ . 
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1. Introduction 

As the title suggests, we are concerned with the recovery of an unknown potential 
of the full line Schrodinger operator — d% + q(x) from a certain set of scattering data. 
In its classical formulation, the scattering data consist of the reflection coefficient 
R, bound states {—Xn}n=i an d their norming constants {c„}^ =1 . 
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2 ODILE BASTILLE, ALEXEI RYBKIN 

It is well-known (see e.g. [9]) that the inverse problem 

S:={R(k),-x 2 n ,c n } kmi ^ N — ► q(x) (1.1) 

for all x is uniquely solvable through the famous Faddeev-Marchenko inverse scat- 
tering procedure. In fact, there is an explicit formula, referred to sometimes as 
Bargmann or Dyson, (see e.g. [10] ) 

q(x) = -25|logdet(l + M a: ), (1.2) 

where M^ is the so-called Marchenko operator, a Hankel integral operator whose 
kernel is constructed in terms of S. 

In practice, however, norming constants {c n } are not available. Moreover the 
authors are unaware of their physical meaning either. Notice that the inverse 
problem 

{R(k)} ke R — > 80*0 (1-3) 

is solved uniquely only if the operator —d 2 + q(x) has no (negative) bound states. 
The classical example of the so-called one-soliton potential 

q(x) = —2m: 2 sech 2 ax + log 

V c / 

produces the scattering data 

{0,-^ 2 ,c} 

with R(k) = for all k, suggesting that the inverse problem (|1 .3[) is ill-posed in 
general. 

The inverse problem (jl.l[) was originally posed and solved for g's decaying suffi- 
ciently fast at infinity. The complete treatment of this problem is done in [2] under 
the Faddeev condition^] 

(l + \x\)\q(x)\dx<oo. (1.4) 



While the Faddeev condition (ll.4[) is typically satisfied in nuclear physics, many 
interesting inverse problems (e.g. in geophysics) deal with potentials which do not 
decay at one of the infinities but still decay at the other infinity. Such potentials 
are commonly called steplike. A suitable analog of inverse problem (jl.ip is well- 
posed as well (see e.g. [1] and the literature cited therein). New circumstances 
arise, however, due to a richer negative spectrum of — d 2 + q(x). But the classical 
Marchenko theory can be successfully adjusted to such setting too. 

In [7] the inverse problem was solved for q(x) — > C, a nonzero constant, as 
x — > — oo and q(x) — ¥ as x — ¥ oo sufficiently fast (with some gaps fixed in [S]). 
The case of g's periodic on the, say, left hand side and decaying on the other have 
been treated first by Hruslov [13] (see also [I]). 

The more general case of g's with no certain pattern of behavior at -co and 
identically zero on the right hand side has been recently treated by one of the 
authors in [TH]. To the best of our knowledge, in the context of steplike potentials, 
the determinant formula (|1.2[) is not available in the literature. 



Certain aspects of the theory developed in [9] actually require the stronger condition 

/ (l + |z| 2 ) \q(x)\dx < oo. 



INVERSE SCATTERING PROCEDURE 3 

The situation with steplike potentials is similar to the case (|1.4j) in that the 
reflection coefficient alone does not determine the potential uniquely. It is natural 
to ask what can make up for the unavailable data in (|1.1|) related to the negative 
spectrum? This problem has generated considerable interest (see e.g. O |4j |6j [12j 
[T5]). In the context of Faddeev potentials, the inverse problem 

{R(k),q(x)} keRx > Q — > q{x),x<0 (1.5) 

is shown to be well-posed. The problem (jl.51) is referred to as the inverse scattering 
problem with partial information on the potential. In |llj , the inverse problem (ll.5|) 
is solved for essentially arbitrary potentials as long as the reflection coefficient can 
be suitably defined. The actual procedure of solving (|1.5|) in [IT] is not scattering 
but rather spectral through solving the Gelfand-Levitan integral equation. 

In [TB] , one of the authors found a way to adapt the classical Marchenko inverse 
scattering algorithm to solve (|1.5|) . No analog of the determinant formula (|1.2p 
appears to be found for inverse problems with partial information on the potential 
(jl.ip in the context of stcplike potentials. The present paper intends to deal with 
this issue. More specifically, we show that if q is locally square integrable, 

( x 2 
sup / |g| < oo, (1.6) 

i<Oii-l 

and the known part q + = g|(o,oo) is subject to (|1.4[) then 

9 (x) = -2<9|logdet(l + (l + M+)- 1 G 2: ) , x < 0. (1.7) 

Here M+ is the Marchenko operator constructed in terms of the scattering data for 
g_l_ and G x is the integral Hankel operator with the kernel expressed in terms of the 
difference R — R + of the (right) reflection coefficients R for the whole potential q 
and R + corresponding to g+ . 

We emphasize that the main issue here is the existence of the determinant in 
(|1.7[) in the classical Fredholm sense. We prove that the operator (1 +M+) G x 
in (|1.7|) is trace class for every x. Our arguments are based upon the simple but 
important observation that the difference R — i?+ is an analytic function (even 
though R and R + are not) and certain limiting procedures (which, as frequently 
happens in similar situations, are quite involved). 

From the geophysical point of view our situation is related to reflection seismol- 
ogy where one is concerned with recovering certain properties q of deeper layers 
of the Earth given already known properties g+ of shallower layers and measured 
reflections R. Of course if the media do not tend to be homogeneous as the depth 
increases, a larger range of frequencies is required to investigate deeper layers. 

Another example of the inverse problem with partially known steplike potential 
is related to neutron reflectometry (see e.g. 2 ) where properties of an unknown 
material are studied by scattering neutrons and measuring their reflection. The 
role of q+ is played by a layer of known material, called a coating. 

The paper is organized as follows. Section 2 lists our notation. Section 3 intro- 
duces the scattering quantities in our setting. Section 4 details relevant properties 
of the Titchmarsh-Weyl m-function which are then applied in our context in Sec- 
tion 5. Section 6 and 7 respectively deal with trace class and Marchenko operators 
and Section 8 contains the main result. 
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2. Notation 



We adhere to standard terminology from analysis, namely R± := [0, ±00), C is 
the complex plane, 

C+ = {z £ C : Imz > 0} , iR+ = {z £ C : z = iy , y £ R+} , 

in the upper half plane, 

R + ih = {z £C : z = x + ih , x £M} 

is the real line shifted h units up. \\-\\ x stands for a norm in a Banach (Hilbert) 
space X. We use (SCI and 5" will typically be R or R±): 

• the usual Lebesgue spaces 

LP(S):=lf:\\f\\ LP{S) :=(^J s \f(x)\ p dx S j P <oo\ , 1 < p < oc 

L°° (S*) := {/ : ||/|| L oo (S) := esssu Pa;e5 |/ (z)| < 00} , 
L Lc ( S ) : = i nLP (A) : A C S 1 , A compact } . 

• the short range or Faddeev class (important in scattering theory) 

L\ (S) = {/ : \\f\\ L{{s) := jT (1 + |*|) |/ (*)| <to < 00 J . 

• the Birman-Solomyak spaces (1 < p < 00) 

r°(Z/(R±)) := 1/ : ||/||^ (LP(R±)) :=jmp ^" _1 |/(x)rd») < 00 1 . 

Next, 62 denotes the Hilbert- Schmidt class of linear operators A: 

e 2 = {A:\\A\\ 62 :=[t r (A*A)f 2 <^} 
while 61 is the trace class 

61 = {A : P|| Si := tr [(A*4) 1/2 j < 00} . 

Spec (A) stands for the spectrum of an operator A and if it is selfadjoint, 
Spec ac (A), Spec d (^4) denote, respectively, the absolutely continuous and discrete 
components of Spec (A). 

The following portion of notation will be used extensively in reference to the 
potential q and quantities associated with it. If xs ( x ) is the characteristic function 
of a set S, i.e. xs (x) — 1, x £ S and otherwise, then we define: 

q+{x) := q(x)xR + , q-(x) = q(x)xR_ , q{x) = q(x)X[- a ,a] fOT some a > 0. 

We also denote 

Sq:=q-q 

and when the cutoff approximation is taken to infinity, i.e. a — > 00, we write Sq — > 0. 
Any quantity X of arbitrary nature (functions, operators, etc.) related to q will 
be denoted X and 

SX := X -X. 
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3. WEYL AND SCATTERING SOLUTIONS OF THE SCHRODINGER EQUATION 

Throughout this section we assume the following. 

Hypothesis 3.1. The potential q is real, locally integrable and such tha\j 

(1) q is limit point case at — oo 

(2) g + is Faddeev class (short range) and hence has the Jost solution at +oo. 

Condition |T]) means that the equation 

— d 2 u + qu = k u (3-1) 

has a unique, up to a multiplicative constant, solution 4 f _, called Weyl, such that 
Sb-(x, k) € L 2 (M_) for any k 2 e C + . Condition ([2]) implies that q is limit point 
case at +oo and the Weyl solution ^+(x, k) can be taken to have the asymptotic 
behavior: 

$+{x,k) = e lkx +o{l) , x->oo (3.2) 

for all reajj k. 

In particular, we have as in classical scattering theory Spec ac (— d 2 + q + ) = 
K_l_. Furthermore, v I / +,5'+ are both solutions to (j3.ll) for a.e. real k and linearly 
independent with constant Wronskian 

W := W U>+{x : k), V+(x, k)\ = 2ik. (3.3) 

So they form a basis of solutions for (|3.1[) for a.e. real k and in particular: 



C(k)V-(x,k) =V + (x,k) + R(k)-$+(x,k) (3.4) 

with some C(k),R(k). We call R the reflection coefficient from the right incident. 
Under our hypothesis, neither C nor R can be analytically continued into the upper 
half plane. Figure [T] illustrates a potential from our hypothesis and the asymptotic 
behavior of C$_ at +oo. 




<1+ = <z|k+ is known 



Figure 1 . Scattering channels for q = g_ + q + 

We now consider separately scattering solutions corresponding to q±. I.e., first 
by our hypothesis at — oo, there is a solution (f-(x, k) to 

:>2„. , „ „. _ ,2, 



-d x u + q_u = k u 



For terminology used without explanation, see e.g. |91|19|. 

the Weyl solution >!/+ coincides in this case with the Jost solution. 
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of the form: (k e K) 



\D(k)*-(x,k) 
\ e - %kx +R_(k)e lkx 



, x < 
, x > 



with some D(k),R-(k) (see FigureEJ. 




FIGURE 2. Scattering channels for q- 



For q+, there exist particular Jost solutions (pi+ and <p r ,+ to 



-c^w + q + u = k 2 u 



such that: (fc e R) 



T+(fc)yv,+ (a;, fc) 



L+(fc). 



-. — ika 



T + (k)^ + (x,k) 
T+{k)e- ikx 



y + (x,k)+R + (k)V+(x,k) 



x < 
a; > 

a; <0 
a; > 



for k real and where T + is called the transmission coefficient and L + , R + the reflec- 
tion coefficients from the left (right) incident. Because q + is short range, T + can 
be analytically continued in the upper half plane and has only a finite number of 
poles {i*e+}n = i- Since, in addition, q + is supported on the right half line, L + can 
also be analytically extended to the upper half plane and shares the same poles as 
T+. However, in general, R+ can not be extended off the real axis. The asymptotic 
behavior of these solutions is illustrated in Figure G3 
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L + {k)e- 

< 



) T + (k)c 



-ikx 



+ R+{k)e lkx + o(l) 



T. 



+(k) e - tkx r 




(1+ 



Figure 3. Scattering channels for q^ 



Note that all right reflection and transmission coefficients can be expressed in 
terms of Wronskians. Of particular interest, we have for a.e. real k: 

W = 2ik = W(^ + (x,k),C(k)^^(x,k)) 



R(k) = 
#_(*) = 


W \V+(x,k),C(k)V 


-{x,k)) 


W 
W(e- lkx ,ip^(x,k)) 
W 


x>0 

o r ^ + (x,k)\ 


R+(k) = 

1 


W (* + {x,k),T + (k)<, 


W 
W(tpr t+ (x, k),lfi. + (x,k)) 



x>0 



T+(k) W 

Note that any truncation q = q\[- a . a ] is compactly supported which implies that 

R, R+ can be analytically continued into C + [9J except at a finite number of poles. 
These poles are located on iR + such that their squares correspond respectively to 
the discrete spectrum of —d\ + q and — c^ + q + . 

4. The Titchmarsh-Weyl to-function 

In this section, we review properties of the Titchmarsh-Weyl m-function which 
is the logarithmic derivative of the Weyl solutions ^±(x, k) as x — > ±0. It will 
be a central object in redefining scattering quantities in the next section. We will 
have to impose some additional conditions on the potential q. Most of the material 
already appeared in [15] but are repeated here for the reader's convenience. 

Definition 4.1. The Titchmarsh-Weyl m-function is defined by: 

d x V±(x,k) 



m±{k z ) = ±- 



*±(s,fc) 



x=±0 



Some of the important properties of the m-function are (see e.g. [14l US]): 

• m± is analytic for all k 2 e C + and has the Herglotz property, i.e. m± : 

C+ -> C+ 

• symmetry m±(z) = m±(z) 
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• the singularities of m± correspond to the spectrum of the half line Dirichlct 
Schrodinger operator, i.e. —d 2 +<jon R± with u(±0) = 0. 

• the Borg-Marchenko uniqueness theorem: m\ = mi => q\ = qi- 

The following representation of m± will be useful. 

Proposition 4.2. Let q be a real function onR such that q e £°°(L 2 (M_))ni 1 (M + ). 

Let 7 = max(7_, 7+) where 

7- = max(y / 2||g_||^ (L2(R _ )) ,e||g_|| foo(L2(R _ )) J , 

Wl+W L 1 (R + ) 

7+ = 2 ■ 

Then for k = a + ih, h > 7, 

»±oo 



/•ztoo 

m±(fc 2 ) =ikT e ±2lkx A±(x)dx 

Jo 



(4.1) 



with some real function A±{x), called the A-amplitude. The integral in (14. ip is 
absolutely convergent and the A-amplitude has the following properties. 

(1) A± — q± is continuous on R± and for ±x > 0: 

\A±(x)-q±(x)\<(±J*\q±(s)\ds} e ±2 ^ (4.2) 

(2) Ifq 1 ,q 2 e£ 00 (L 2 (R-))nL 1 (R + ) then 

q x {x) = q 2 (x) on [0,±o] =4- A x (x) = A 2 (x) on [0,±a], (4.3) 

(3) for any h > j, 

||e=F 2te ^ ± (a ; )|| il(R±) <C(/ l ,g ± )<oo 

and C(h,q±) is a nonincreasing function of h. 

(4) For any n > 7, 

e 2ha: v4_(x)eZ. 2 (R_). 



Proof. The representation (14. ip appeared in [14] for short range g's and in [18] for 
more general g's. Properties (I4.2I) - (|4.3I) were derived for q + £ L 1 (M.+ ) in [T5] then 
for g+ £ £ 00 ( J L 1 (K + )) in [5] but since 

m_(g_(x),fc ) = m + (g_(— x), A; ) 

and £°°(.L 2 (R_)) C ^ 00 (£ 1 (R_)), we have adjusted the results accordingly. So only 
([3}-(j4]) require a proof. We will consider A- and p — 1,2. Using Minkowski's in- 
equality, one needs only to show e 2hx q-(x) and e 2hx (A-(x) -g_(x)) are in £ P (R_). 
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Dropping the subscripts, we have 



pO 

\e 2hx q{x)\ P dx = / e 2hpx \q{x)\ p dx 





oo ._„ 


e 2hpx \q(x)\ p dx 


< 


■y^ e -2hpn 

oo 


f " \q(x)\ p dx 

J-n-l 


< 


V^ e -2hpn 
n=0 


\\q{x) £oo( L p( R _)) 




1 


?W u»rr,pm m 



For the next term, we will make use of the following: (x < 0) 

V x \ /•- n+l ~N 

" l#»(H(]R_)) 



/■0 -H ,-n+l -N 

/ |«(*)|<fe< £ / \q(s)\ds<J2h\\ 

J x „ i J — n „ i 



< (1 #) IMI.«°°(.L 1 (K-)) — ^ X ^ IMI^°o(£p(R_)) ■ 

where the last inequality is a direct consequence of Holder's inequality. So, 

,o 

\e 2hx {A{x)-q{x))\ P dx= / e 2hpx \A(x) - q{x)\ p dx 

J — OG 

2p 



< /" ftoOx-i)* ( J \q{s)\ds\ dx 

< T e ^-7)x (1 _ a;) 2p|| g ||2p r ^ 



|£°°(Lp(R_)) 
^ — oo 

One readily verifies that for any m = 0, 1, 2, . . . and b > 

11' 



-oo " ,,_ n 



Therefore, we have 

3=0 

So ©-gl) are verified with 



loo 1 ^ V ^ ^ "' - [2«(/,- 7 )l 2p+1 ^ 



■21:., „.., ,..,,!'- .,..^ 2 ^ ! Y^ [2p(ft-7)]" || ,|2p 



K^-J ~ J _ e -2h H9-WIU~(i 1 (R-))" , "J7^I ) 3 ^ j! l^-H^CiHR-)) • 

Similarly for A + , we can take in ([3]) 

1 2 

C(h,q+) = \\q+\\ L i {R+) + -^—^ lk+llz,i( K+ ) ■ Q 

Remark 4.3. In the case of a truncated potential q, since 7± < j± and C{h, q±) < 
C(h,q±), all above results remain true for the same h. If, in addition, q + € 
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L L( R +) then 9+ e t°°(L 2 (R + )) and thus e - 2hx A + {x) G L 2 (R+) for h large 
enough. 

Corollary 4.4. Let q € £°° (L 2 (K-)) , anrf Zei /i > 7- where 7_ is defined as in 
Proposition \4-%\ Then 

(1) ik-m-{k 2 ) <= L 2 (R + ih), 

(2) i/gL = g_|r_ 0) oi /or some a > 0, tten 

6m-(k 2 ) := m_(fc 2 ) - fh-(k 2 ) ->• m L 2 (R + i/i) , a -> 00. 
Proof. Note that for fc = a + i/i where ael,A>7_, 



ifc — m_(fc ) 



e- Aai V' ,x A_(z)da; 



where by Proposition 14.21 e 2hx A-(x) e L 2 (R + i/i). The Plancherel formula in our 
case takes the form: 



-2ia:c 



f(x)dx 



= V^\\f(x)\\ L 2 {R ^ 



LH 



and hence 



Ufc — 771_(fc ) 



L 2 (R+i/i) 



= V^ \\e 2,lx A.{x) 



L 2 (R_) 



< OO 



In the case g_ = g-|[- a ,o] ; we have 7- < 7- and by Proposition ^. 21 (|2j). 



<5m_(fc ) 



e- Mkx 5A(x)dx 



and therefore for any h > 7_ , 

||5m_(fe 2 )|| r2 ,„, ..s = V7r||e 2/la:! (Jvl(x)|| r2 ,, n -* , a ->• 00. 

II v ' II L z (R+ih) v II v ' II L 2 {{ — 00, — aj) 



D 



5. The reflection and transmission coefficients 
In this section we establish some properties of one of our main objects: 

G(k) := AJZ(ft) = R{k) - R+(k). 

As mentioned when first introduced, neither R nor R + can be analytically ex- 
tended to the upper half plane for a potential q under the conditions of Hypothesis 
13.11 (or those in Proposition 14.21) . But by rewriting G exclusively in terms of R- , 
T+, and L+, we will see that G can be analytically extended to the upper half 
plane. We also derive key properties of R- , T + , L + and G in C + which will be used 
later to recover q-(x) - assuming R, R + are known. 

First we rewrite the reflection and transmission coefficients in terms of the m- 
function. Setting 



rn± = m±(k + iO) := lim m±(k +ie), 

e->0+ 
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wc have: 



R(k) = 
R+(k) = 
R-(k) = 
L + (k) = 
T+(k) = 



m_ 


+ 


m + 


*+(o,ft) 


rh~+ 


+ 
+ 


m+ * + (0,fc) 
ik$ + (0,k) 



ik — Ta- 



lk + m_ 
ik — m+ 
ik + m + 



2ik 



(5.1) 
(5.2) 
(5.3) 
(5.4) 
(5.5) 



(jfc + TO + )* + (0,fc)' 

The above are obtained using continuity of the various solutions and their deriva- 
tives in x at the point x = or, alternately, the Wronskians. Recall that the above 
are defined for a.e. real k. 

From (|5.3p - (|5.4|) and properties of the Titchmarsh-Weyl m-function, one notes 
that L+, R- have a meromorphic extension to the upper half plane. By Proposition 
4.21 |3]), L + , R_ are smooth on R + ih for any h > 7. Furthermore, by the Borg- 
Marchenko uniqueness theorem L + , i?_ determine uniquely respectively q± (x) [14j . 

Now, using (|3.3|l and (|5.1|l - (|5.5|l . one readily verifies that: 



G(k) 



(5.6) 



l-L+(k)R-(k) 

and thus G can be analytically extended to the upper half plane (recall that T. 
also has an analytic extension since q + is short range [S]). 

Proposition 5.1. Let q be as in Proposition ^.'^ and let 

h± = inf < h : h > -f± and C(h, q±) < — 

Then for all h > max(/i + , hJ) 

(1) R-,L+ £ (L°° n L 2 ) (R + ih) with their L°° norm no greater than 1/3. 

(2) 8L+ -t in L°°(R + ih) when a -> 00. 

(3) R- G L^R + i/i). 

(4) 5R- -> mi^R + i/i). 

Proof. The above results are direct consequences of Proposition 14.21 For all k £ 
R + ih, such that ft. > max(/i+, /&_), 

■>±oo 



ife — m±(k ) = 



i/j + m ± (/c 2 )| = 



(■: 



±2ifex 



:x At(x)dx 



< 


\e* 2hx A ± (x)\\ L 


ft. 

i(K±) - 2 




pzkoo 


= 


2ikT / e 2lkx A ± (x)dx 


> 2\k\ ■ 


1 

1 2|fc| 


/ e ±2 > 


> 


3|fc| 







c ,4 ± (x)dx 
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2< 



where we have used |fc| > h. Thus jfc+m t k i\ € L (K + ift) and it follows immedi- 
ately from dSU-dUl} that ||^-(fc)|| L =o (R+ ^) ,\\L + (k)\\ Lao{m+zh) < | and L + ,i?_ G 

L 2 (R-M/i). 

We further obtain i?_ G L 1 (K + i/i) by the Cauchy-Schwartz inequality for 
h > h- : 



\R- 



L 1 (R+i/i) 



< 



1 



ik + rri-(k 2 ) 



\ik — m_(fc ) 



L 2 (M+ih 



L 2 (R+ih) 



By Remark 14.31 the above is also true for R_, L+. 
Now for all k G R + i/i with h > h + , 

—2ik5m + (k 2 ) 



\5L + (k)\ 



(ik + m + (k 2 ))(ik + m + (k 2 )) 



<- h \6m +{ ki)\ 



<- « znx 8A + (x) L,,. ,, -> , a^oo. 



We also have 



SR-(k) = — 



-2ik5m,-(k 2 



(ik + rri-(k 2 ))(ik + m_(fc 2 )) 
Using I/°° norms and the Cauchy-Schwartz inequality: 



r-^-lli^R+i/i) - 



-2ik 



4 
< - 
~ 3 



ik + m,-(k 2 ) 
1 



L°°(R+ih 



<5m_(fc 2 ) 



i/c + m_(& 2 ) 



L 1 (M+ih) 



ifc + m-(k 2 ) 



L 2 (R+ih) 



■ ||*"»-(* 2 )|L»(R+ 



!7l) 



and the right hand side, by Corollary 14.41 goes to zero when a — > oo. 



□ 



Remark 5.2. Property ([3]) um7Z J^ay a crucial role. Note that if q(x) — c5(x), 
where 6 is Dirac 's 8-function, then 

2ik — c 

which is not in L 1 (M. + ih). This suggests that the condition £°°(L 2 (R)) may not be 
relaxed to read £°°(I 1 (R)). 



Corollary 5.3. For any finite z , andq-,h under the conditions of Provosition [5A 

e 2lkz R-(k) G L 1 (M. + ih). 



Proof. Immediately follows from e R-(k) \\ Tl 



L 1 (R+ift) 



e- 2 '-\\R_(k)\\ Ll(M+lh) . 



D 



While trivial, the above corollary plays an important part in our arguments. The 
reflection coefficient for the shifted potential q(x + z) is R(k)e 2lkz where R(k) is 
the reflection coefficient corresponding to q(x). 

Lemma 5.4. Let q + G L\(M.+) and let h > (3 where 

/3 = 2maxj||<7 + || L i (R+) ,lj. 
Then 



\T + (k)\\ 



L°°(R+ih) ' 



T+(k) 



L°°(R+ih) 



<2? 



and ST + (k) -> in L°°(R + ih) as a -*■ 00. 
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Proof. The following are well-known facts (e.g. [9]) for q+ short range and sup- 
ported on R+ : 

(1) T + (k) is analytic in <C + except at a finite number of simple poles {iXn}n=i 
where 



N < 1 



\x\ \q+(x)\dx. 



(2) |T + (fc)T + |L+(fc)| = 1 for a.e. real k and T+(-k) = T+(k). 

(3) T + (k) admits the following representation for any k £ C + : 



r + (*)-n^«pfi 

We also have the Lieb-Thirring inequality [3j 

N 



logiT+nr 1 



id — k 



duj 



X^rt ^ L i/2.i / |«+(a;)|da; 

i ili 



where 1/2 < L 1/2 ,i < 1.005. Thus, for any jfe G K + ih, h > j3 



\T + {k)\ 



N 

n 

n=l 



fc + i>d 



f£ I^Cfi 



exp | — / Re ■ 



■log|T+(w)| X doj 



n^+ 



4/i-+ 



^ expf^/.^'y^U 



Since |T + (w)| < 1 for a.e. real w, the above becomes 

jv / 



\r + (k)\ < n ji 



4ftxiJ 



71=1 



(h - xtf 



But by (|5.8p . we also have for each n: 



^n < ^1/2,1 lk+ll L i(R+) < o/ 3 < n h 



0/2 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



and hence for all k £ R + ih where h > j3 

\T + (k)\ < (±±) < 2? 

where we have used N < (3 from (|5.7p . Inequalities (|5.7[) - (|5.10[) are also valid for 
T + (k) and thus T+,T+ are uniformly bounded on M + i/i. 
From [SJ, we now use the following results 

1 If 

7F77T = 1 - ^rr / g+(o;)j/ + (x, fc)dx 

r+(fc) 2ifcy R 

where y+(x,k) := e~ (fi t +(x,k) satisfies for Imfc > 0, 



y+(x,k) = l+ / D k (t-x)q+(t)y + (t,k)dt, D k (y) = 
|y + (:r,fc)l<^(/3)(l + M) 



2ifc : 



(5.11) 
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where if is a constant depending only on j3. Note that 
T + (k)T+(k) 



ST +^ 2lk 

From (J5TT) . 



q+(x)6y+(x,k)dx + / y+(x,k)5q+(x)dx 



(5.12) 



5y+(x,k)= / D k (t-x)y + (t,k)5q + (t)dt + D k (t- x)q + (t)5y + (t,k)dt (5.13) 

J X J X 

and since for k ^ 0, |Dfc(2/)| < m - for all y > 0, Imfc > 0, 



K(ft) u I, ri«+(*)i 

^II^ + IIlj( R+) +/ -^ 



|*l/+(x,A)| < ^T Z II^ + II L ;«^+ / ^F 1 !*»+(*, fc)|«« 



i^(/3) fl 

<^rre^||^+|| L j fE , , fc^O , Lmfc>0 



t)y_l_ (x,k)k 

Hence for f|5.12p . we have 



by iteration on the Volterra integral equation for k /q^ s i i - — derived from (J5.13D . 



\\ST + (k)\\ L ^ R+lh) <2^^1\\5q + \\ L{iw+) 

where the right hand side goes to zero for a — > oo by the dominated convergence 
theorem. □ 

Proposition 5.5. Under Hvvothesis \3.1[ 

AR(k) := R(k) - R+(k) 
is analytic in C + except on a set 

S = {ix+}{J<rCM+ 
where 

{- K) 2 } = Spec d (-^ + q+) , a = {\:\ 2 eSpec(-d 2 x + q)nR-}. 
Proof. By direct computation, we have 

m + (k z ) + m^(k z ) 
where 



l + L+(k) * + (0,*) 
We gather the following facts: 

(1) it is well-known that T + (k),L + (k) are analytic in C + \ {i*tn} n =i wnere 

< — (x+) > is the negative simple discrete spectrum of —d 2 + q+. So we 

also have that g(k) is meromorphic in C + but with poles different from 
those of T + , L + . Poles of g(k) correspond to the poles of m+, i.e. n's such 
that *+(0,k) =0. 

(2) recall that m± is analytic in C + except for some singularities k G ?R+, 
hence so is (^ "+^" vfc 2 ) ' Note that m + (k 2 ) + m_(fc 2 ) = corresponds 



Singularities of m_|_ are a finite number poles since q_|_ is short range I14| whereas the set of 
singularities of m_, while bounded by 7_, need not be made of poles - it could be continuous. 
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to W (ip~, ip+) = 0. But then if these two solutions to <9 2 + q are linearly 
dependent, then ip±(x, k) £ L 2 (R) and so A G a. 
(3) we now consider singularities of m±. Those cases correspond exactly to 
ip±(0,K) — with ?/>j_(0,k) 7^ 0. But if t/>_(0, k) = 0, we can assume 
iI>+(0,k) 7^ (otherwise k £ a) and so to + (k 2 ) finite. Therefore, by (|5.14p . 

T(k) 

G(K) = -i^) 

is finite unless k £ {ix^}. Now if ip + (0,K,) = 0, we have by (J5.5I) that 

rr ( \ 2lK 

T + (k) 



is finite and by (|5.14[) since we can assume to_(k 2 ) is finite, then 

G(K) = T+W v;(o,«) 

is finite too. 
Thus, we find that G(fc) is analytic in C + \ ({i)t+} U cr). D 

Remark 5.6. 5y Proposition \4-ty we have 7 > sup |<S| and thus G{k) is smooth 
onM. + ih for any h > 7. -/Vote £/ia£ cr need not be finite, but is bounded. 

Proposition 5.7. Let q be a real function on R smc/i i/ioi g G ^°°(L 2 (]R_))nL}(R + ) 
and Zei h > h where ho = max(n + , /i_ , /3) where h± are as in Proposition \5.1\ and 
f3 as in Lemma \5.4\ Then 

G(k)=AR(k)£L lr R + ih) and SG(k) -> m L 1 (R + i/i). 



Proof. Omitting the variable fc for brevity, from f|5 . 6[) . Proposition 15 . 1 1 and Lemma 
I5.4[ we have: 

\\G\\ L ^ + ih) < f ■ 2 2/l ° • WR-hip+ih) < °°- 
By direct computation, 

SG=Gl 6L + + G 2 6T + + G 3 5R- 
(I - L+R_)(l - L + R_) 
where 

Gi = T 2 i?_fi_ , G 2 = R-{1-L + R_)(T + + T + ) , G 3 = ? 2 . 
By Proposition 15. 1[ 

1 



(1 - L+RJ){1 - L + RJ) 



< 

L°°(i+a) 



so it is enough to show that each term inside the brackets in (I5.15[) goes to zero in 
L 1 (R + ih) as a — > 00. Indeed 

2 2ft, 

T 



l^l^+W ^(R+ih) - ~T~ ' W-^-WL^R+ih) ' II<^+IIl°°(B+j/i) ~> 0' 



\\G 2 ST + \\ L1(m+ih) < — ■ 2 ° ||i?-|| L i( K+ j ft ) • \\^ T +\\ L ^(R+ih) "*■ °> 

|G3<5i?-|| L1(R+lh) < 2 2,l « • ||^-|| L1(R+ ^ ->• 0. D 
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Corollary 5.8. Let z be a fixed real parameter and let q, h be as in Proposition 
\5l\ Then G z {k) ■- e 2tkz AR(k) € L 1 {R + ih) and 5G z {k) ^ in L l (R + ih). 

Proof. Note that G z (k) correspond to the shifted potential q(x + z). For such 
potential, R-(k) becomes R-(k)e 2lkz , T+{k) remains the same and L+{k) becomes 
L(k)e~ 2tkz . So by Corollary [O] and Proposition EH G z (k) £ L l (R+ih) and 
6G z (k) -> in L^R + ih). □ 

6. A TRACE CLASS OPERATOR 

In this section we introduce a lemma which appeared in a more general form 
in |17) and will be a central argument in the main result of this paper in the next 
section. 

Proposition 6.1. Let A be smooth onM. + ih for some h > 0, and A £ L 1 (M. + ih). 
Then the integral operator A on L 2 (M + ) with kernel 

dk 



&{x,y) 



Jk(x+y) 



A{k) 



K+ih 



2vr 



x,y > 



is trace class, and 



IIAIk^lHl^cW 



Proof Denote Ah(a) = A(a + ih) and f{z) = ^ J M e lkz f{k)dk. Then rewrite A as 
an operator on L 2 (M) by considering x, y £ R and: 

A(x,y) = x(x)e- Hx+y) A h (x + y) X (y) 

where \ is the characteristic function on R + . By convolution and a change of 
variable, we have: 



So 



A h (x + y) = ( y/A h * sJA h j(x + y) 

/~Ah{x - s)\f~Ah{y + s)ds. 
A1A2 where Ai, A2 are operators on L 2 (R) with kernels 
Ai(x, s) = X {x)e~ hx yfM{x - s), 
A a (*,y) = x{y)z~ hy ^A~h{y + a). 



One readily has 



62 



\A k ((;,n)\ 2 d£d V 



X(z)e 



1 
4nh 



- 2hz dz 

2 



>A h 



1 



dS 



LH 



Anh 



\A h \ 



where we have used the Plancherel equality 



2^11/H 2 and hence 



< 



1 



le 2 ll A 2|l S2 - |^ II^IIlI(K) - 4^ W A \\Li(R +l h) 



D 
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7. Classical Marchenko inverse scattering 

In this section we review well-known facts about how to recover a potential q in 
the Faddeev class from the scattering data associated with the Schrodinger operator 
—d x + q(x) via the Marchenko inverse scattering procedure (see e.g. [Tll9lll0|). 

For q € L\(M), the scattering data consisting of 

• the discrete spectrum {— x^}n=4 °f the Schrodinger operator —d x + q{x) 



onL 2 ( 

normi 

operator, 



• norming constants {c n }^ =1 associated to the bound states of the Schrodinger 



• and the reflection coefficient R(k), k <G R 
determine together the potential uniquely. By the inverse scattering procedure 

q(x) = -2d x k x (0 + ), (7.1) 

where k x <E L 2 (R + ) solves the Marchenko equation 



k x (y) + M x (y) + M x (y + z)k x {z)dz = , y>0 (7.2) 



with 



M x {-)=M(- + 2x), (7.3) 

A i r 



M(s) = ^ cie~^ s + — J e^R{k)dk. (7.4) 



71=1 

If we define the Marchenko operator on L 2 (R + ) as: 

(M x f)(y)= M x (y + s)f(s)ds , f e L 2 (R + ) (7.5) 

Jo 

then (|7.2p becomes 

(1 + M x )k x {y) = -M x (y) 
and 1 + M x is boundedly invertible [9] . 

Assuming the Fredholm determinant in (|7.5|) is well-defined, one can also rewrite 
dLH as [10] 

q(x) = -2dl\ogdet(l + M x ), (7.6) 

known as the Bargmann, Dyson, or determinant formula (see e.g. |16j). The 
determinant is well-defined if M^ is trace class. However, we don't know if it is the 
case for a generic short range potential. 

We choose to detour this fact. To this end, we express the Marchenko kernel 
in a different form. Recall the well-known fact (see e.g. p]) that if a short range 
potential is supported on K_, then R(k) can be analytically continued in C + , its 
poles are {i>c n }^ =1 and 

Res(R(k),iK n ) = %c 2 n , 

where c„ is the norming constant associated to n n in the scattering data. 

So for any h > max{xr n }, one can deform the contour |16j in (|7.4[) and by the 
residue theorem rewrite the Marchenko kernel (17.41) 



M{s) = tt [ e iks R{k)dk. 



(7.7) 



Note that (I7.7P can then be used for any compactly supported q's with compact 
support using a shifting argument. 
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Remark 7.1. We will consider a potential which is locally square integrable on 
the line, in £°° (I? (R-)) and such that q+ is Faddeev class. The classical inverse 
scattering results do not apply directly since q is not short range and the negative 
part of the spectrum of —d 2 + q(x) need not be finite. However, since q e L^ oc (R) c 
Lj oc (M), we have 

• q G £i(R) so the classical inverse scattering procedure applies to the trun- 
cated potential, 

• q is compactly supported so we can use (|7.7[) , 

• q G ^°°(L 2 (R)) which we will show implies that M. x is trace class and so 
(|7.6p applies. 

The above will be a basis for our limiting procedure. 

8. Main result 

We now present our main result which gives a formula to recover a nondecaying 
unknown potential g_ assuming that g+ and the reflection coefficient R are known. 

Theorem 8.1. Let q be a real, locally square integrable potential on R such that 

(q± = q\m±) 



• sup / |g_(s)| ds < oo (uniformly in Lf oc ) 

x<OJx-l 

• I (1 + x) \q+(x)\ dx < oo (short range) 



and let R{k), R+(k) be the right reflection coefficient corresponding to q,q+ respec- 
tively. 

Let M+ be the Marchenko operator associated with the scattering data 

[R + (k),-{x n ) >Cn} fceR)1 < n <jv 

for q+ (given by (|7.3|) - (|7.5|) ) and let G x be the Hankel integral operator associated 
with R — R + . I.e. 

(P*f)(v)= G x (y + s)f( S )ds , feL 2 (R+), (8.1) 

Jo 

G x (s) = ^- f e lk(s+2x \R~R + )(k)dk (8.2) 

Z7T JR+ih 

with some h > sufficiently large. 
Then for any x < 

q_(x) = -2d 2 log dct (l + (1 + M+) -1 G x ) (8.3) 

with the determinant defined in the classical Fredholm sense. 

Remark 8.2. Theorem \8.1\ solves the inverse scattering problem for a steplike 
potential with the knowledge of its short range part. Indeed, given (short range) q+ 
one solves the direct scattering problem and finds the scattering data 

{R+{k),-{*n) 2 ,C n } km ^ n < N 

forq+. Then, we construct M+ by (|7.3[) - (|7.5p and given the (right) reflection coef- 
ficient for the whole potential q, one constructs by (|8.ip - (|8.2p the Hankel operator 
G x . The unknown (non decaying) part of q is recovered for each x < by (|8.3|) . 
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Proof. We will first prove the statement for q. For a fixed a > 0, q is compactly 
supported. Hence, R can be analytically continued in C + except at a finite number 
of poles {i>(n}n=iJ an d the Marchenko kernel (|7.4I) becomes 



M(s) = i 

Z7T 



e R(k)dk , ft, > max{xr„} 



A' 



+ih 



Dehne q a (x) 
£°°(L 2 (R_)). 



<j(x + a) as in Figure HI Then q a is supported on M_ and g Q e 




Figure 4. Shifted potential g a (^) = q(x + a) 

Now its reflection coefficient i? a ._(fc) G ^(R + ih) for h > h a where /i a is defined 
as in Proposition 15.11 for q a . But R(k) = R a .-(k)e~ 2tka , so they share the same 
poles and by Corollary [OJ R(k), R{k)e 2lkx e L 1 (R + ih) for h > h a . We also have 



that R(k)e 2lkx is smooth onl + i/i for h > max{x„}™ =1 so by Proposition ^. 1[ M x 
is trace class. Hence the following Bargmann formula applies: 

q{x) = -25|logdet(l + M a: ). (8.4) 

Now write 

R = R+ + G , G = Ai? = R — ii_j_ 

and split the Marchenko operator accordingljO 

M x =M+ + G X . 

The same /i a is enough to ensure M+ e ©i since R+(k)e 2l,ka corresponds to q a (x) 
(above the shaded region in Figure [3| . 

In addition, for h > ho where ho is the samqj as in Proposition 15.71 G,G,SG € 
L 1 (R+i/i) and by Proposition l5.5l and the subsequent remark, we also have G, G, 5G 
smooth onl + ih. Since by Corollary |5.8[ the same applies to G x , G X ,5G X , we can 
apply Proposition 16. II and conclude that G x , <G X , SG X £ 6i. 

Therefore, first we rewrite the Bargmann formula (|8.4[) as: 

q(x) = -2dl logdet(l + M+ + G x ) 

= -23£logdet(l + M+)-23glogdet(l + (l + M+) 



where we have used the fact from classical Marchenko theory that 1 
edly invertible. But 

q+(x)=-2d 2 logdet(l+M+) 



X G X ) (8.5) 

+ M T isbound- 



5 Because q+ is compactly supported, M+ can be equivalently expressed as 117.41 1 or 17.71 1. 
Note that ho is independent of a. 
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and q+{x) — for x < 0. So (|8.5|) becomes for x < 0: 

g_(a:) =-2^1ogdet(l + (l + M+)- 1 G K ). (8.6) 

Now, we now use the fact that (1 +M Z ) remains boundedly invertible in its limit 
so the right hand side of (|8.3p is well-defined. In addition, by Proposition l6.ll for 
h > h 

\\8 < &x\\ 6l ^ ^^ \\?> G x\\ L-LQiL+ih) 

and the right hand hand side of the inequality goes to zero by Corollary 15.81 for 
a — > oo. Therefore, we find indeed that the limit of (|8.6I) is (|8.3|) . D 



Note that if M+ G ©i then (153)) simplifies to 

q(x) = -2<9jlogdet(l + M !C ), iel, (8.7) 

where M x = M+ + G x . It is, of course, well-known (see e.g. [9]) that under our 
condition on g+, M T G ©2 but we couldn't prove it for ©1. We were unable to 
find a rigorous proof of such a statement in the literature either. (It is typically 
assumed (frequently implicitly) or referred to as "too involved".) However, since 
M x G © 2 , then dct(l + M x ) can, in fact, be regularized differently from (|8.3[) (sec 
Q7] for details). 

In conclusion, we emphasize that the fact that (|8.3[) is understood in the classical 
sense is indeed quite important as it guarantees the convergence of various types 
of approximation of (1 + M+) G x in trace norm. This, in turn, means a certain 
stability of the inverse problem algorithm based upon (|8.3p . 
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